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Abstract 



A manifestly gauge invariant formulation of chiral theories with fermions 

on the lattice is developed. It combines SLAC lattice derivative Q, [P, Q and 

\ generalized Pauli-Villars regularization Q . The theory is free of fermion dou- 

■ bling, requires only local gauge invariant counterterms and produces correct 

CN . results when applied to exactly solvable two dimensional models. 

(N 

a^ 1 Introduction 

In this paper I discuss a possibility to describe chiral fermions on the lattice in 
the framework of the SLAC formulation HH, 0, supplemented by the additional 
^ , Pauli-Villars (PV) type regularization which allows to avoid the problems of the 

original SLAC model. 

k>( ; It is known that due to nonlocality of the SLAC action perturbation theory suffers 

^ • from the appearance of nonlocal counterterms and spurious infrared divergencies . 



The attempts to overcome this difficulties by making a partial resummation when 
applied to solvable models lead to the wrong physical spectrum 0|. 

The origin of these difficulties is a singular behavior of the fermion propagator 
and vertices near the edge of the Brillouin zone ^ ^- The contribution of this 
region produces nonlocal counterterms and infrared divergencies. Essentially the 
same problem is present in all other models proposed to eliminate fermion doubling 
on the lattice p|, [|], ||10|- (For detailed references see [III], ||l3l)- In particular 



in Wilson formulation the contribution of the region near the edge of the Brillouin 
zone produces gauge noninvariant counterterms. In other words a lattice formula- 
tion is not a bona fide regularization for chiral fermions, and to get an invariant 
calculational scheme one needs to suppress a contribution of momenta of the order 
of cut-off. In fact it is the same problem which has to be solved if one wishes to 
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construct a gauge invariant regularization for the chiral fermions in the continuum 
theory. 

Recently we proposed a manifestly invariant regularizaion for the anomaly free 
chiral models, in particular for the standard model |^ (see also IjT^). Our procedure 
is a generalization of a Pauli-Villars regularization which allows for an infinite num- 
ber of auxiliary PV fields. It has been shown that when applied to the chiral models 
on the lattice this procedure leads to a gauge invariant continuum theory without 
fermion doubling both in the case of Wilson fermions JTSf and in the Smit-Swift 
model However both these formulations have some drawbacks. It is a lack of 
a gauge invariance for a finite lattice spacing in the first case, and the necessity to 
introduce an additional Yukawa interaction in the second case. It seems that the 
most appropriate way to implement this regulaization in lattice models is to use the 
SLAG formulation. 

We shall show that if one writes a lattice regularization for anomaly free chiral 
gauge models following |@|, @ and introduces simultaneously PV fields according 
to [0,the resulting model is manifestly gauge invariant, has no fermion doubling, 
and does not require nonlocal counterterms in perturbation theory. When applied to 
solvable two dimensional models this procedure leads to a correct continuum limit. 



2 Four dimensional SO{10) model. 

In this section we consider the grand unified 50(10) model. It is worthwhile to 
emphasize that representations we consider are not real, and the choice of 5*0(10) 
is explained by the possibility to write all the formulae in a compact form. We 
could consider equally well the standard model or Weinberg-Salam model provided 
the lepton and quark representations are choosen in such a way that anomalies are 
compensated. The SLAG lattice action for the grand unified 50(10) model looks 
as follows: 

1= E ^X{x)^,iD,{x~y)Pe^v{i E 9M^)}^liy) (1) 

k,fi,x,y Zfj,=Xn 

Here = A^^aij. The matrices aij are the 50(10) generators:crjj = l/2[ri,rj], 
where Fj are Hermitian 32 x 32 matrices which satisfy the Glifford algebra: [Fj, F^] = 
25ij. The chiral 50(10) spinors ip± = 1/2(1 ± Fn)?/), where Fn = F1F2 . . . Fio, 
describe the 16-dimensional irreducible representation of 50(10) including quark 
and lepton fields. We assume also that the spinors ip± are left-handed ip± = 1/2(1 + 
75)V^±- Index k numerates different generations. D^{x) is the SLAG derivative 

= (^^^/^ exp{iA;x} (2) 

The action (|l]) is manifestly invariant with respect to 50(10) gauge transformations 
and does not suffer from fermion doubling. However it is not local which leads to 
difficulties we mentioned above. 
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To avoid these problems one modifies the action (|I]) by adding the analogous 
interaction of PV fields 0], [0, [|T§ 

I^I + Ipv (3) 

Mr - 

Ipv= 'ilJrix)-ff,iD^{x - y)Pexp{i J2 9At,iz)}i>riy) - ^Aix)CDCTnipJ+ 

M - - 

+ ^rix)-f^tTiiD^{x-y)Pexp{t ^ gA^{z)}^riy) - ^Mx)CdC^J + h.c. 

Here ipr are anticommuting PV fields and $r are commuting PV fields, which real- 
ize the reducible 32 dimensional representation of 5*0(10), Cd is the usual charge 
conjugation matrix and C is the 50(10) conjugation matrix aJ^C = —Caij. 

The action (|^) is manifestly invariant with respect to the 50(10) gauge trans- 
formations. Now we shall show that it generates a perturbation theory which in the 
limit a ^ results in a continuum theory with only local gauge invariant countert- 
erms. Nonlocal counterterms do not appear. 

The perturbative propagators look as follows 

P 

S{[,+^+ = -p^ (4) 

P 

r, _r, _r,_ _r. _ MrCpCTu 



P2 + M2 



where P^ is the sawtooth function 



P^{p) = - 2m-- (2m- 1)- < (2m + l)-. (7) 
a a a 

The expansion of the action (^ in terms of g generates the interaction vertices with 
the increasing number of vector lines. (For detailes see f^, Q). The three point 
vertex looks as follows 



with 



= expl^oA^ 
ia 

One sees that the interaction vertex (^ is nonlocal, and it introduces additional 
singularities into Feynman diagrams. Indeed, if p < ^, but g > ^, the interaction 
vertex becomes ~ ^ which results in the appearance of the nonlocal divergent 
contributions 0. 

When the PV fields are introduced the contribution of the region p ~ is 
suppressed and nonlocal divergencies do not appear. 
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Let us consider for example the vacuum polarisation diagrams generated by the 
vertex (^). They look as follows: 

n?/^-^/!|;,Tr[Mi±r,,Mx (10) 

Tr [ (1 + 75)7 iPil) + M )7 [Pa + p)+M . Ml)-PAl + P)]m+p)-m] 

(We consider the contribution of the chirality preserving propagators (^ ^)). The 
contribution of the original fermion is given by n+i°), {Mr = 0). Contribution of the 
bosonic fields differs by sign. 
First of all we note that 

Tr [aijTuCTki] = 0. (11) 

and therefore positive and negative " chirality" P V fields give the same contribution 
to H^u- Secondly, the nonlocal contribution comes from the region of integration 
where - — \Pfj,\ < < -■ li one chooses -C a~^, one can expand the integrand 
in this region in terms of Mj.. Terms of zero order produce the total contribution 
proportional to 

A; + 2^(-l)V, (12) 

r 

Here Cr is the number of the PV fields with the mass Mr. The contributions of 
fermionic and bosonic fields differ by sign. The factor 2 arises because there are PV 
fields of both chiralities, and the factor k is due to the presence of k generations of 
original fermions. If the number of generations is even one can always choose the 
coefficients in such a way that 

fc + 2^(-l)''c, = (13) 

r 

Imposing the further PV conditions 

J2cr{-irM^ = (14) 

r 

one can make the contribution of the domain ^ — < |/^| < ^ vanishing in the 
limit a — * 0. In the remaining integral over |/^| < ^ — 

P^{l)-P^{l+P)r~.p^ (15) 

Therefore the nonlocal factors K^{p) are compensated and in the continuum limit 
one gets a manifestly gauge invariant expression for the polarization operator. 

The same arguments are obviously applicable to the diagrams including chirality 
changing propagators (||). 

In the case of an odd number of generations the eq.(pr3|) cannot be satisfied for 
any finite number of PV fields, as by construction the coefficients Cr are integer. As 
was discussed in ref. M, |]T3[, in this case the desired suppression of the region 



near the edge of the Brillouin zone may be achieved if one introduces an infinite 



4 



series of PV fields with tlie masses = Mr, M{a) —>■ oo when a 0, M <^ a ^. 
In this case after summation over r one gets the integrand 

, R' = P^{l) + P^{l) + Pm + P3il) (16) 



MRsinh{nRM~^) 

which vanishes exponentially for i? ~ ^. Therefore this region does not contribute 
and one can use the same arguments as above to show that in the continuum limit 
one gets a manifestly gauge invariant polarization operator and only usual local 
gauge invariant counterterms are needed. 

The proof for the diagrams with 3 and 4 external lines is given in the same way. 
The diagrams generated by the vertex function and the propagators (|D,(||) may 
be presented in the form 

T ■^^ll...^ln (^^\ 

where 

Ij^^:,. = -g"- Jl ^ Tr ± TuWi • • • x (18) 

Tr [ (1 + 75)7^1 {P{1) + Mr)i,, . . . {P{1 - Pn) + Mr)] X 

[P,,{1) - P,Al + Pl)] . . . [P,Al-Pn) - P,M 
[P\l)+M^]...[P\l-Pr,)+M^] 

Due to the properties of V matrices Tr [ai^jj^n . . . <Ji„j„] = for n < 4 and as before 
the contributions of positive and negative chirality fermions are equal. The integrals 
over dli may be separated into two parts: 



dli + / dli, Vin = < \pi\, Vout = \k\ > \Pi\ (19) 

Vin J Vout (X 

Here \pi\ = sup \p} + . . . pl\, I = 1, . . . n — 1. If belongs to V^n then the factors 
Pf^{l + p^ + . . . p'-) — P^{1 + p^ + . . . p'~^) can be replaced by p|^,and the integral over 
Vin reduces in the limit a — > to the usual gauge invariant continuum expression. 
Hence if we adjust the parameters of the PV fields in such a way that the integral 
over Vout vanishes in the continuum limit, our statement will be proven. 

In complete analogy with the analysis of the polarization operator in the case 
of even number of generations of the original fields we may take a finite number of 
PV fields with the parameters satisfying eqs. ( ]TB| , |Tip. Expanding the integrands in 



(|T^) one sees that the leading terms cancel and the remaining terms vanish in the 
continuum limit. 

In the case of an odd number of generations one needs an infinite system of PV 
fields. Taking the traces and separating tensorial structures in eq. (|18D one can 
write for 1^^...^,^ 

~ fl dH +2? n^i Af-^'"(/,Q,M,) 
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where Af^'"'^" is a polinomial in M^. The summation over r can be done exphcitly 
as in ref. 1^, |T^. One gets 



p dH if-""(/,g,) 

Mi...Mn ~ ^27r)4 Mi?sinh(7ri?M-i) ^ ^ 

In the region Vout the integrand vanishes exponentially when a — 0, and the integral 
over Vin reduces to the gauge invariant continuum expression. 

To extend this proof to the higher order diagrams we firstly note that contrary 
to the statements which one can meet in the literature the model (|l]) generates 
a finite number of types of divergent diagrams. Indeed, scaling the integration 
variables — >• l^a~^, in eq.( |l^) one would get for arbitrary diagram the highest 
degree of divergency ~ a~^. However this estimation is too rude. Let us show 
that all the diagrams with n > 5 are convergent. To prove it we separate the 
domain of integration overdli as in eq.( |T^). If belongs to Vm then the factors 
P^j(/ + + . . . p'') — P^{1 + p^ + . . - p^'^) can be replaced by pjj. Each such factor 
decreases the degree of divergency by one. The integral over Vout is limited by 
the product of independent external momenta pip2 ■ ■ -Pn-i- Indeed, the integrand 
contains the factors 

Pi{l + api + ... + apk) - Pi{l + api + . . . + apk+i) (22) 

This factor can be different of apk+i only on the interval of values of which is 
equal to ap^+i. Noting that the remaining terms in the integrand are not singular 
at p\j^i = 0, we conclude that the integral is limited by F\apl._^_-^\ where F is a 
polynomial over ap^+i- It follows that the diagram with n external lines contributes 
the factor 

< a"1ap^il • • • \ap^^„-l\ (23) 

and therefore all the diagrams with n > 5 are convergent. 

These reasonings are easily extended to the case when the diagram includes the 
higher order vertices r„(/,pi . . .p„). The vertices r„ are defined up to nonessential 
factors by the recurrent relations 

Tn+l{l,Pl ■ ■ - Pn+l ~ -r„(/,pi ...Pn)+ Tn{l + Pn+l,Pl ■ ■ ■ Pn) (24) 

Let us consider some momentum p*^.It follows from eqs.( ^, |2^) that is the sum of 
Pi depending on different arguments and all the terms in r„ may be separated into 
pairs whose arguments differ by p^. Therefore r„ is proportional to p^ everywhere 
except for the interval of the values of U whose length is equal to ap^. In complete 
analogy with the case discussed above we conclude that the eq. (^) holds. Counting 
the number of independent momenta one sees that the diagrams with n > 4 are 
convergent. In the limit a — > they coincide with the gauge invariant continuum 
expression. The corresponding diagrams with PV fields vanish in the continuum 
limit. 

It is worthwhile to notice that these reasonings do not exclude nonlocal singu- 
larities for diagrams with n < 4. For example in polarization operator there is only 
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one independent momentum and without PV regularization it contains nonlocal 
singularities. 

This discussion shows also that spurious infrared divergencies which could be 
caused by the singularities of the vertex functions are in fact absent being supressed 
according to the eg. (|23|) . 

The analysis of the diagrams with divergent subgraphs in particular overlapping 
divergencies requires more complicated machinery and has not yet been done. We 
just mention that this problem may be avoided by introducing an additional higher 
derivative regularization for Yang-Mills field, as it has been done in the refs. |r5|, 



Tql . With this regularization all the above results are trivially extended to arbitrary 
diagrams. 

It completes the proof of the absence of nonlocal counterterms in our scheme. 
Of course all the arguments given above used the weak coupling expansion and it 
would be very important to verify them nonperturbatively. Having this in mind we 
consider in the next section two dimensional chiral models. 



3 Two dimensional models 

In two dimensional models the procedure described above simplifies greatly because 
the only divergent diagram is the second order polarization operator. In the contin- 
uum Abelian model the polarization operator is the only nonzero diagram and for 
that reason the anomaly free chiral models are essentially equivalent to the vectorial 
Schwinger model. It is not true in the nonabelian case. 

We start with the anomaly free chiral Schwinger model. The fermionic part of 
the action is 

^ = i'^+i^h^i^D^ix - y) exp{iea ^ A^{z)}^p+{y) + (25) 



^!(x)7^iL'^(x -?/)exp{z— =a 

k=l,2 v2 Zf,=Xf, 



where 



/■f (fk 

D/^ix) = / " -——ikf,exp{ikx} (26) 
(27r)2 

-0-1- = |(1 ± 75)?/^. The coupling constants of the positive and negative chirality 
fermions are adjusted in such a way that the 75 contributions cancel and the model 
is anomaly free. 

To suppress the contribution of the region near the edges of the Brillouin zone 
we introduce the vectorial interaction of bosonic PV fields $ 

I = J2[^i^hf^^D^{x-y)exp{tea Af,{z)}(!>{y) + M^{x)^x)] (27) 
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The regularized action = I + Ipy is invariant with respect to the gauge transfor- 
mations 



tp+{x) ^ eic^{ieai{x)}i)+{x)] iptix) ^ e^^{i-^ai{x)}^-{x) (28) 

<l>(x) ^ exp{iea{(a;)}$(x), A^{x) A^{x) + a'^f^x^ + a) - f{x)]. 

Consider firstly the vacuum polarization graph. The corresponding amplitudes 
for each fermions are 

n£^" - /! Tr [ 1(1 ± 75)7.P(07.P(/ + V)] x (29) 

[P^{l)-P^{l+p)][P^{l+p)-P^{l)] 
P^{l)P^{l + p) 

(30) 

Here ± stands for the contributions of positive and negative chirality fermions , and 
(7+ = e, (7- = The corresponding amplitudes for the PV fields are 

= Jl Tr [7,(P(/) + Mrh^{P{l +p) + Mr)] X (31) 

[P^(l) - P^{1 + p)][P^(l + p) - P^{1)] 

[P2(/)+M2][P2(/+p)+M2] • 

The complete contribution of the original fermions is given by 

n;:+ + 2np + n;^+ + 2nJ; (33) 

The terms proportional to 75 cancel and the sum ( |33D coincides up to the sign and 
mass terms with the PV amplitudes (^, ^). Following the procedure described 
in the preceeding section we can separate the integration domain into Vm and Vout- 
In the domain Vout, p ^ M and one can expand the integrands in terms of M. 
The first terms cancel and the next terms are majorated by M'^a? . Therefore the 
integral over Vout vanishes in the continuum limit. The integral over Vin in the limit 
a — >■ coincides with the PV regularized continuum theory. The same reasoning are 
applied when Iq belongs to Vin and h belongs to Vout and visa versa. 

The diagrams with more than two photon lines are analysed in the same way 
as above. One can easily see that the infrared divergencies are absent. Indeed, 
according to the estimate given above the integrals over fermion loops are limited 
by the product of independent external momenta. Therefore the infrared singular 
factors are supressed. A straightforward power counting shows that the higher order 
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sea-gull vertices give a contribution vanishing in the continuum limit. For example 
the tadpole diagram produced by the second order vertex contributes the factor 



I ^ K,{k)K,{k) ^-^^^ 

where -D^^ is the photon propagator. This expression is infrared singular for p ~ 
-. However it was shown in the preceeding section that the integration over p 
produces the factors ~ A;^ compensating this singularity. Therefore one can rescale 
the integration variables fc^ a~^k^, which allows to estimate this integral when 
a ^ 0: 

n ~ a / £xF{x, ap) (35) 



where it is understood that 11 is attached to some fermion loop with the integration 
momentum p. According to the discussion above the integral is convergent and 
therefore 11 — when a — 0. Analogous reasonings show that all loops including 
internal photon lines vanish in the continuum limit. The only nonzero diagram is 
the photon polarization operator which coincides with the continuum result: 

n,. = --(V-^) (36) 

7r 



A;2 



Note that if one would apply to our model Rabin procedure without additional PV 
regularization one would get for the mass gap the result which is the continuum 
one times \/2 [0]. To avoid misanderstanding we emphasize that we did not solve 
exactly the Schwinger model for a finite lattice spacing, but rather demonstrated 
that the lattice approximation developed above generates in the continuum limit 
the correct exact solution of the Schwinger model if the limit is taken termwise in 
the weak coupling expansion. 



4 Discussion 

In this paper we showed that modification of the SLAC model by introducing the 
generalized PV regularization provides a manifestly gauge invariant lattice descrip- 
tion of anomaly free chiral models. Our model does not suffer from fermion doubling 
and leads in the continuum limit to the usual gauge invariant results with only local 
counterterms. When applied to the anomaly free chiral Schwinger model it repro- 
duces in the continuum limit the well known exact solution. 
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